变角度纤维复合材料层合斜板的颤振分析

假设纤维方向角沿层合板的长度方向线性变化，研究了变角度纤维复合材料层合斜板的颤振．通过坐标变换将斜板变换为正方形板，采用层合板表面连续变化的速度环量来模拟空气对其的作用，速度环量分布利用Cauchy积分公式计算．建立了系统的Lagrange方程并采用Ritz法得到了层合板的自振频率和颤振/不稳定性分离临界速度．通过数值算例验证了本文模型和方法的正确性和收敛性，分析了各个铺层内纤维方向角的变化对自振频率和颤振/不稳定性分离临界速度的影响．研究结果表明，通过纤维的变角度铺设，可有效地提高层合板的基频和颤振/不稳定性分离临界速度．经合理设计的变角度复合材料层合板具有抑制颤振的作用．     

机翼大攻角下失速颤振的气动弹性研究

本文采用流固耦合的数值方法研究了机翼在0°～50°攻角下的颤振。计算结果表明,随着来流攻角α0的增大,机翼的固有频率对颤振的影响越来越大,颤振由线性的强迫振动逐渐发展成为非线性的自激振动,而当α0增加到一定程度以后,大尺度分离流交替地从机翼头部和尾部产生,机翼和流场会发生共振,引起机翼的失速颤振。     

用“参数多项式方法”确定风力透平叶片的失速振荡流

用“参数多项式方法”确定风力透平叶片的失速振荡流陈佐一,孙永忠,杨玲（清华大学热能工程系北京１０００８４）关键词：参数多项式方法,风力透平,颤振在现代各类流动诱导振动问题的研究中,流体激振的安全性分析与诊断,要求在较短的时间内,对大量的复杂流动工况进...     

The chaotic response of a fluttering panel: The influence of maneuvering

The influence of maneuvering on the chaotic response of a fluttering buckled plate on an aircraft has been studied. The governing equations, derived using Lagrangian mechanics, include geometric non-linearities associated with the occurrence of tensile stresses, as well as coupling between the angular velocity of the maneuver and the elastic degrees of freedom. Numerical simulation for periodic and chaotic responses are conducted in order to analyze the influence of the pull-up maneuver on the dynamic behavior of the panel. Long-time histories phase-plane plots, and power spectra of the responses are presented. As the maneuver (load factor) increases, the system exhibits complicated dynamic behavior including a direct and inverse cascade of subharmonic bifurcations, intermittency, and chaos. Beside these classical routes of transition from a periodic state to chaos, our calculations suggest amplitude modulation as a possible new mode of transition to chaos. Consequently this research contributes to the understanding of the mechanisms through which the transition between periodic and strange attractors occurs in, dissipative mechanical systems. In the case of a prescribed time dependent maneuver, a remarkable transition between the different types of limit cycles is presented.Nomenclature a plate length - a _r u _r /h - D plate bending stiffness - E modulus of elasticity - g acceleration due to gravity - h plate thickness - j₁,j₂,j₃ base vectors of the body frame of reference - K spring constant - M Mach number - n 1 + ₀/g - N ₁ applied in-plane force - p–p aerodynamic pressure - P pa ⁴/Dh - q ₀/2 - Q _r generalized Lagrangian forces - R rotation matrix - R ₄ N, a ²/D - t time - kinetic energy - u plate deflection - u displacement of the structure - u _r modal amplitude - v₀ velocity - x coordinates in the inertial frame of reference - z coordinates in the body frame of reference - Ka/(Ka+Eh) - - elastic energy - 2qa ³/D - a/_mh - Poisson's ratio - material coordinates - air density - _m plate density - - _r prescribed functions - _r sin(r z/a) - angular velocity - a/v₀ - skew-symmetric matrix form of the angular velocity     

结构非线性颤振半主动抑制

本文利用数字仿真技术，对结构非线性颤振半主动抑制－颤振驯化方案进行了探讨，仿真结果表明，对于一定的结构非线性类型和参数，利用非线性颤振的极限环特性，可使系统的颤振响应被抑制在幅值很小的稳定域内，从而达到减缓颤振的目的。     

Fluid Flow-Induced Nonlinear Vibration of Suspended Cables

The nonlinear interaction of the first two in-plane modes of a suspended cable with a moving fluid along the plane of the cable is studied. The governing equations of motion for two-mode interaction are derived on the basis of a general continuum model. The interaction causes the modal differential equations of the cable to be non-self-adjoint. As the flow speed increases above a certain critical value, the cable experiences oscillatory motion similar to the flutter of aeroelastic structures. A co-ordinate transformation in terms of the transverse and stretching motions of the cable is introduced to reduce the two nonlinearly coupled differential equations into a linear ordinary differential equation governing the stretching motion, and a strongly nonlinear differential equation for the transverse motion. For small values of the gravity-to-stiffness ratio the dynamics of the cable is examined using a two-time-scale approach. Numerical integration of the modal equations shows that the cable experiences stretching oscillations only when the flow speed exceeds a certain level. Above this level both stretching and transverse motions take place. The influences of system parameters such as gravity-to-stiffness ratio and density ratio on the response characteristics are also reported.     

Chaotic Analysis of Nonlinear Viscoelastic Panel Flutter in Supersonic Flow

In this paper chaotic behavior of nonlinear viscoelastic panels in asupersonic flow is investigated. The governing equations, based on vonKàarmàn's large deflection theory of isotropic flat plates, areconsidered with viscoelastic structural damping of Kelvin's modelincluded. Quasi-steady aerodynamic panel loadings are determined usingpiston theory. The effect of constant axial loading in the panel middlesurface and static pressure differential have also been included in thegoverning equation. The panel nonlinear partial differential equation istransformed into a set of nonlinear ordinary differential equationsthrough a Galerkin approach. The resulting system of equations is solvedthrough the fourth and fifth-order Runge–Kutta–Fehlberg (RKF-45)integration method. Static (divergence) and Hopf (flutter) bifurcationboundaries are presented for various levels of viscoelastic structuraldamping. Despite the deterministic nature of the system of equations,the dynamic panel response can become random-like. Chaotic analysis isperformed using several conventional criteria. Results are indicative ofthe important influence of structural damping on the domain of chaoticregion.     

机翼/外挂颤振主动抑制的控制律研究

本文介绍了对机翼/外挂系统的颤振抑制控制律的研究。并用三种控制律的结果进行了风洞实验验证。结果表明:控制律的设计是成功的,理论计算与实验结果吻合良好。 文中重点介绍了两种控制律,它们都是以现代控制理论为基础。首先,把最优控制理论与颤振分析的状态空间法相结合,得到状态反馈。然后导出不同的输出反馈。文中,还讨论了该系统的阵风减缓和稳定裕度。 为对比起见,还给出了在同一模型上,用气动能量概念方法导出的控制律。由此,可看出它们在颤振抑制效益上的差异。     

Deteriorated Geometrical Stiffness for Higher Order Finite Elements with Application to Panel Flutter

Higher order elements were first design for linear problems where, in certain situations, they present advantages over the lower order elements. A method to efficiently extend their use to geometrical nonlinear problems as panel flutter and postbuckling behavior is presented. The chaotic and limit-cycle oscillations of an isotropic plate are obtained based on direct integration of the discretized equation of motion. The plate is modeled using the von Karman theory and the geometrical nonlinearities are separated in a nonlinear term of the first kind which manifests especially in the prebuckling and buckling regimes, and a nonlinear term of the second kind which is responsible for the postbuckling behavior. A fifth order, fully compatible element has been used to model thin plates while the inplane loads where introduced through a membrane element. The aerodynamics was modeled using the first order 'piston theory. The method introduces the concept of a deteriorated form of the second geometric matrix which is equivalent to neglecting higher order terms in the strain energy of the plate. This allows for a drastic reduction in the computational effort with no observable loss of accuracy. Well established results in the literature are used to validate the method.     

Experimental investigations on the dynamic behavior of a 2-DOF airfoil in the transitional Re number regime based on digital-image correlation measurements

The present paper investigates the fluid–structure interaction (FSI) of a wing with two degrees of freedom (DOF), i.e., pitch and heave, in the transitional Reynolds number regime. This 2-DOF setup marks a classic configuration in aeroelasticity to demonstrate flutter stability of wings. In the past, mainly analytic approaches have been developed to investigate this challenging problem under simplifying assumptions such as potential flow. Although the classical theory offers satisfying results for certain cases, modern numerical simulations based on fully coupled approaches, which are more generally applicable and powerful, are still rarely found. Thus, the aim of this paper is to provide appropriate experimental reference data for well-defined configurations under clear operating conditions. In a follow-up contribution these will be used to demonstrate the capability of modern simulation techniques to capture instantaneous physical phenomena such as flutter. The measurements in a wind tunnel are carried out based on digital-image correlation (DIC). The investigated setup consists of a straight wing using a symmetric NACA 0012 airfoil. For the experiments the model is mounted into a frame by means of bending and torsional springs imitating the elastic behavior of the wing. Three different configurations of the wing possessing a fixed elastic axis are considered. For this purpose, the center of gravity is shifted along the chord line of the airfoil influencing the flutter stability of the setup. Still air free-oscillation tests are used to determine characteristic properties of the unloaded system (e.g. mass moment of inertia and damping ratios) for one (pitch or heave) and two degrees (pitch and heave) of freedom. The investigations on the coupled 2-DOF system in the wind tunnel are performed in an overall chord Reynolds number range of $9.66\times 10^3\le \text{Re}\le 8.77\times 10^4$. The effect of the fluid-load induced damping is studied for the three configurations. Furthermore, the cases of limit-cycle oscillation (LCO) as well as diverging flutter motion of the wing are characterized in detail. In addition to the DIC measurements, hot-film measurements of the wake flow for the rigid and the oscillating airfoil are presented in order to distinguish effects originating from the flow and the structure.